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Abstract 

Searches for first generation scalar and vector leptoquarks, and for squarks in 
R-parity violating SUSY models with the direct decay of the squark into Stan¬ 
dard Model particles, have been performed using e’''e“ collisions collected with 
the OPAL detector at LEP at e+e“ centre-of-mass energies between 189 and 
209 GeV. No excess of events is found over the expectation from Standard Model 
background processes. Limits are computed on the leptoquark couplings for dif¬ 
ferent values of the branching ratio to electron-quark hnal states. 
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1 Introduction 


The observed symmetry between the lepton and quark sectors in the Standard 
Model is not yet understood, but could be interpreted as a hint for common under¬ 
lying structures. Consequently, many extensions of the Standard Model postulate 
the existence of leptoquarks (LQ), which are coloured spin 0 or spin 1 particles 
carrying both baryon (B) and lepton (L) quantum numbers. The Buchmiiller- 
Riickl-Wyler (BRW) model [|I| adopted in this paper assumes lepton and baryon 
number conservation. Two additional assumptions on the leptoquark couplings 
are made in the following: only the couplings A within one generation of fermions 
are assumed to be non-zero (to respect lepton flavour conservation), and only the 
case of chiral couplings is considered, i.e. it is assumed that the product ArAl 
of couplings to left-handed and right-handed leptons vanishes (to respect lepton 
universality). With the latter assumption the branching ratio /3e of leptoquarks 
to final states with a charged lepton is restricted, as shown in Table |I|. 

This paper presents the results of a search for single production of leptoquarks 
in e^e“ collisions with the OPAL detector at LEP. Leptoquarks would be pro¬ 
duced by an elect ron-quarkQ fusion. In the dominant diagram, which is shown 
schematically in Figure |I], the quark is the result of a process where one of the 
incoming leptons radiates a photon which subsequently fluctuates into a hadronic 
state 0. Under the above assumptions on the leptoquark couplings, only first 
generation leptoquarks coupled to electrons could be produced in electron-photon 
scattering, and they would decay into either an electron and a quark or into a 
neutrino and a quark. Therefore the main signature of single leptoquark events is 
one hadronic jet balanced in the transverse plane either by one isolated electron 
or by missing transverse energy due to the neutrino. The hadronic photon rem¬ 
nant would disappear down the beam-pipe or add some activity in the forward 
region of the detector. 

In supersymmetric models with R-parity violation, scalar quarks (squarks) 
have the same production mechanism as some leptoquarks. The topology of 
events with single squark production depends on the squark decays; if R-parity 
conserving decays can be neglected, the results on the left-handed couplings for 
leptoquarks presented in this paper apply also to the corresponding squark states. 
In the following, leptoquarks and the corresponding squark states with R-parity 
violating decays are generically referred to as leptoquarks. 

Experiments at LEP |^, |^, HERA |^, and the Tevatron have searched 
for leptoquarks. We present a search for leptoquarks with Mlq > 80 GeV in 

^Charge conjugation is implied throughout this paper for all particles, e.g. positrons are also 
referred to as electrons. 
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scalar LQ (q) 

charge 

F 

coupling and 
decay mode 

/3e 

So (or dfl) 

-1/3 

2 

Al : e^u, i/Ld 
Ar : ej^u 

1/2 

1 

So 

-4/3 

2 

Ar: e^d 

1 

Si/2 (or cIl) 

+1/3 

0 

Al: 

0 

Si/2 (or ul) 

-2/3 

0 

Al: ehd 

1 


+2/3 


Al : i^lU 

0 

Si 

-1/3 

2 

Al : r'Ld, e^u 

1/2 


-4/3 


Al : e^d 

1 




Al : i^LU 

0 

Sl/2 

-2/3 

0 

Ar : ej^d 

1 


-5/3 


Al : ClU 

1 




Ar : ej^u 

1 





coupling and 


vector LQ 

charge 

F 

decay mode 

/ 3 e 


-1/3 


Al : i^Lcl 

0 





Vl/2 

-4/3 

2 

Ar : ej^u 

Al : ehd 

1 

1 



Ar: e^d 

1 

Vl/2 

+2/3 

-1/3 

2 

Al : i^lU 

Al : ehu 

0 

1 

Vo 

-2/3 

0 

Al : ehd, r'LU 
Ar : epj^d 

1/2 

1 


+ 1/3 


Al : 

0 

Vi 

-2/3 

0 

Al : ehd, i^lU 

1/2 


-5/3 


Al: e^u 

1 

Vo 

-5/3 

0 

Ar: ej^u 

1 


Table 1: The first generation scalar (S) leptoquarks/squarks and vector (V) lepto- 
quarks in the BRW model according to the nomenclature in with their electric 
charge in units of e and fermion number F = L + 3B. For each possible non-zero 
coupling A the decay modes and the corresponding branching ratio (3e for the 
decay into an electron and a quark are also listed. The restrictions on the values 
of f3e arise from the assumption of chiral couplings. 
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electron-photon scattering using data corresponding to an integrated luminos¬ 
ity of 612.3 pb“^ at e’''e“ centre-of-mass energies between 189 and 209 GeV. The 
data sample is split into eight bins with average e^e“ centre-of-mass energies of 
188.6 GeV, 191.6 GeV, 195.5 GeV, 199.5 GeV, 201.6 GeV, 204.9 GeV, 206.5 GeV, 
and 208.0 GeV. This paper extends the previous OPAL analysis to the high¬ 
est LEP centre-of-mass energies and uses rehned experimental techniques. The 
data used in [Q are included in this analysis and thus the results from @] are 
superseded. 


2 The OPAL Detector 


The OPAL detector is described in detail elsewhere @]. It is a multipurpose ap¬ 
paratus with almost complete solid angle coverage. The central detector consists 
of two layers of silicon micro-strip detectors [^j and a system of gas-hlled tracking 
chambers in a 0.435 T solenoidal magnetic held which is parallel to the beam axis. 


A lead-glass electromagnetic calorimeter with a presampler is located out¬ 
side the magnet coil. In combination with the forward calorimeters, the forward 
scintillating tile counter [^, and the silicon-tungsten luminometer |^, a geo¬ 
metrical acceptance is provided down to 25 mrad from the beam direction. The 
silicon-tungsten luminometer measures the integrated luminosity using small- 
angle Bhabha scattering events il^ . The magnet return yoke is instrumented for 


hadron calorimetry, and is surrounded by several layers of muon chambers. 


3 Monte Carlo Simulation 


The Monte Garlo simulation of the process e^e“ —LQ-t-X-fe and the calculation 
of the production cross-section, dehned as the sum of the cross-sections for the 
particle and antiparticle states, are performed with the program ERATO-LQ [^. 
The hadronisation of the leptoquark decay products and the hadronic photon 
remnant is performed using JETSET [Q. Details of the Monte Garlo simulation 
of single leptoquark production and the calculation of the total cross-section can 
be found in 0. Samples of 3000 signal events are generated for both scalar 
and vector leptoquark states, for both eq and uq hnal states, for e+e“ centre- 
of-mass energies ^/s^ of 200 and 208 GeV, and for scaled leptoquark masses 
xlq = ■'^lq/'v/s^ of xlq = 0.42, 0.53, 0.63, 0.74, 0.85, 0.90, 0.95, and 0.98, 
making a total of 64 samples. For a given scaled leptoquark mass the event 
properties depend only weakly on the centre-of-mass energy. 
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To study Standard Model background processes, two-fermion hadronic events 
(e’''e“ —>• qq(7)) are simulated with PYTHIA 5.722 while the e’''e“ —>• r+r“ 
and e^e“ —*• e^e“ processes are generated with KORALZ 4.02 [^] and BHWIDE [IT 


respectively. Deep inelastic hadronic two-photon background events in the range 
> 4.5 GeV^, including charged current deep inelastic scattering (CC DIS) 
events, are simulated with HERWIG 5.8 [|^, while PRO JET 1.10 is used to 
generate hadronic two-photon events in the range < 4.5 GeV^. For leptonic 
two-photon events, Vermaseren [113 is used. Final states with four-fermion pro¬ 


duction are simulated with grc4f pO 


All Monte Garlo events are passed through a full detector simulation and 
the same reconstruction algorithms as the real data. 


4 Event Analysis 


We search for events with one hadronic jet and either an electron or missing energy 
balancing the transverse momentum of this jet. The analysis uses tracks recon¬ 
structed in the central tracking devices and clusters measured in the calorimeters 
and forward detectors. The selection of tracks and clusters is similar to that 
used in previous OPAL analyses. In addition to quality requirements [||, tracks 
must have more than 20 hits in the central jet chamber, more than half the num¬ 
ber of hits expected given the polar angle of the track and the geometry of the 
jet chamber, and a transverse momentum with respect to the beam direction of 
more than 120 MeV. Galorimeter clusters must pass energy threshold cuts to sup¬ 
press noise. To avoid double counting of particle energies, a matching algorithm 
between tracks and clusters is applied ||22[ . 


Jets are reconstructed from the tracks and clusters using the Durham 
and cone algorithms. The y parameter value yi 2 ( 2 / 23 ) where the transition 
from a one- to two-jet (two- to three-jet) event occurs is determined with the 
Durham algorithm. These parameters provide information on the event shape 
in the following analysis. A cone jet hnder |^, requiring a minimum energy of 
15 GeV in a cone of half-angle 1.0 radian is used to provide information on jet 
directions and energies. No requirement on the number of jets is made. 


4.1 The Electron Plus Hadronic Jet Channel 


Tracks are identihed as electrons if more than 20 measurement samples can be 
used to determine the specihc ionisation energy loss, dE/dx, if the dE/dx prob¬ 
ability |2^ for the electron hypothesis exceeds 1%, and if the ratio of the electron 
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energy measured in the electromagnetic calorimeter to the track momentum is 
between 0.7 and 2.0. In each event, the electron with the largest momentum is 
assumed to be the one from leptoquark decay. 

Candidate leptoquark events are then selected based on the following cuts: 

EQl) The event must contain more than Eve tracks, and the energy measured in 
the hadron calorimeter must exceed 1 GeV. 

EQ2) The ratio .^t/ E^is of the missing transverse energy and the visible energy 
Eyis has to be smaller than 0.2. This cut mainly reduces background from 
4-fermion and two-photon events. 

EQ3) The event must contain an identihed electron, and its energy has to be 
larger than 0.2yG^, where ^/s^ denotes the e’''e“ centre-of-mass energy. 
The energy of the cone jet containing the most energetic electron must not 
exceed the electron energy by more than 12% of the centre-of-mass energy. 
This ensures that isolated, energetic electrons are selected. 

EQ4) At least one cone jet not containing the electron is required, and the lep¬ 
toquark is reconstructed from the electron and the most energetic other 
cone jet (quark jet). The quantity cos6** is the cosine of the helicity angle 
between the electron in the leptoquark rest frame and the leptoquark flight 
direction in the laboratory frame. Since high-mass leptoquarks will be al¬ 
most at rest in the laboratory, the finite energy and momentum resolution 
of the detector will produce an apparent leptoquark direction close to the 
jet or the electron direction. Therefore the distribution of cos^* will be 
strongly peaked at ±1. To remove small-angle Bhabba events in which one 
electron produces a shower in the tracking chambers, the multiplicity of 
the quark jet (dehned as the number of tracks plus unassociated clusters) 
is required to be larger than 10 for | cos6**1 > 0.95 and | cos6*e| > 0.8, where 
6*e denotes the electron polar angle with respect to the outgoing electron 
beamQ. 

The above cuts are applied for all centre-of-mass energies and leptoquark masses 
and states. To reject background further, an XLQ-dependent likelihood is con¬ 
structed, where the leptoquark mass is reconstructed as the invariant mass of the 
electron and the quark jet. The likelihood uses as inputs: 


• the charge-signed cosine of the electron polar angle, —geCos6*e, 

^The OPAL coordinate system is defined as a right-handed Cartesian coordinate system, 
with the X axis pointing in the plane of the LEP collider towards the centre of the ring, and 
the z axis in the direction of the outgoing beam electrons. 
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• the quantity cos 6 ** defined above, 

• the logarithm of j /23 from the Durham jet finder, 

• the cosine of the polar angle of a third jet in the event, if one exists^ 
multiplied by the sign of cos^e, and 

• the ratio ^tZ-EvIs- 

Monte Carlo samples that are statistically independent of those used in the effi¬ 
ciency determination are used to obtain the probability density to construct the 
likelihood function. Two-fermion events, four-fermion events, and two-photon 
events are considered as separate background hypotheses. The signal reference 
distributions for the likelihood for a given event should ideally be derived from 
Monte Carlo simulation of a leptoquark with a mass corresponding to the mea¬ 
sured electron-jet mass of that event. Since the scaled electron-jet mass xlq of 
the event does not in general correspond to one of the signal Monte Carlo masses, 
bin-by-bin interpolation from the two Monte Carlo samples with a:LQ bracketing 
that of the measured event is used to determine the signal reference histograms 
used for its likelihood. For each event, the signal likelihood is calculated from 
the average of the scalar and vector leptoquark distributions at the reconstructed 
leptoquark mass. Finally, events are accepted if they pass a cut on the likelihood 
value C > 0 . 8 . 

In Table ||, the number of data events and the expected number of Standard 
Model background events are shown after each selection cut. The selection effi¬ 
ciencies for two different scalar and vector leptoquark masses at 200 GeV e'''e“ 
centre-of-mass energy are also given. The likelihood interpolation technique de¬ 
scribed above was checked at each leptoquark mass simulated in the Monte Carlo 
by comparing the efficiency obtained when using reference histograms derived di¬ 
rectly from the Monte Carlo at that mass to that obtained when the histograms 
were derived by interpolating the next two nearest simulated mass points. It was 
found that the interpolation method degrades the efficiency by about 15% for 
low-mass leptoquarks, while it has a very small effect for high-mass leptoquark. 
The efficiency for leptoquark detection is therefore reduced to account for this 
effect. 

Figures ^-e show the distributions of some of the cut variables for data. 
Standard Model background and the leptoquark state S 1/2 with two different 
scaled masses xlq. In general, good agreement between the data and Standard 
Model background is observed. The distribution of reconstructed leptoquark 
masses is shown in Figure ^f for events that pass all cuts. No significant excess 
is observed in the data. 

^Events in which no third cone jet is found are collected in a separate bin. 
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selection steps (cumulative cuts) 

events 

(EQl) 1 (EQ2) 1 (EQ3) | (EQ4) | £ > 0.8 


signal efficiencies: 


scalar LQ, ttilq = 106.0 GeV 

72.6% 

55.0% 

39.8% 

39.7% 

(15.6 ± 1.1)% 

scalar LQ, tulq = 190.5 GeV 

79.4% 

69.3% 

57.3% 

56.7% 

(49.2 ± 2.7)% 

vector LQ, tulq = 106.0 GeV 

66.5% 

46.8% 

32.4% 

32.4% 

(9.8 ±0.8)% 

vector LQ, ttilq = 190.5 GeV 

70.6% 

61.9% 

51.5% 

51.1% 

(47.7 ±2.7)% 

expected background events: 

e’''e“ —>■ 2 fermions 

46391 

17682 

94.0 

91.5 

9.1 ± 2.9 

e’''e“ —>■ 4 fermions 

12800 

5111 

431 

428 

23.7 ± 7.4 

two-photon 

765676 

5702 

71.0 

67.3 

11.9 ± 3.9 

charged-current DIS 

25.8 

0.1 

< 0.1 

< 0.1 

< 0.1 

total expected background 

824892 

28496 

596 

587 

44.7 ± 14.0 

observer 

. events: 

data 

735325 

28888 

562 

536 

43 


Table 2: Selection efficiencies in the eq channel. The listed efficiencies are for 
^/s^ = 200 GeV. The remaining number of data events and the expected number 
of background events are also listed after each selection cut. All event numbers 
are quoted after a loose event preselection in which more than 2 tracks, more 
than 0.1 GeV total energy in the hadron calorimeter, and an identified electron 
are required. The errors are the total systematic uncertainties including the 
Monte Garlo statistical error. The component which is not due to Monte Garlo 
statistics is correlated between the separate contributions to the background. The 
determination of the systematic errors is discussed in Section 


10 








































4.2 The Neutrino Plus Hadronic Jet Channel 


In this case we search for events with a single hadronic jet whose transverse energy 
is balanced by missing transverse energy. The kinematic properties of leptoquark 
events decaying to z/q final states do not depend strongly on the leptoquark 
mass, so the selection is independent of the reconstructed transverse mass and 
no likelihood is used. The selection cuts are: 

NQl) As in the eq channel, the event must contain more than five tracks, and the 
energy measured in the hadron calorimeter must exceed 1 GeV. 

NQ2) The presence of at least one cone jet is required. To ensure that the event 
is well contained in the detector, the polar angle of the missing momen¬ 
tum vector must satisfy | cos6*(^| < 0.95, and the polar angle of the most 
energetic cone jet | cos^jetil < 0.85. 

NQ3) The transverse mass of the leptoquark is calculated from the transverse 
momentum of the most energetic cone jet in the event and the missing 

transverse momentum of the event, — cos a), where 

a is the opening angle between and For one-jet events, there is a 
strong correlation between and and the ratio is required to satisfy 
0.497 < < 0.512. The cut is chosen to be asymmetric to retain 

good efficiency for low-mass leptoquarks. 

NQ4) Finally, the requirement of lnp 23 + 1-4 ■ lnpi 2 < —10.5 is used to select 
events with a one-jet topology, as shown in Figure The factor of 1.4 has 
been found to give good separation of signal and background. 


Table || shows the numbers of data and expected background events nor¬ 
malised to the data luminosity after each cut. The efficiencies for two different 
scalar and vector leptoquark masses are also given. 

Figures ^-c show some of the selection variables for the data, the Standard 
Model background and the leptoquark state S 1/2 with two different scaled masses 
xi^Q. The expected background describes the data well. The transverse mass 
distribution of the selected leptoquark candidates is shown in Figure P after all 
cuts. No significant excess is observed in the data. Therefore, limits are calculated 
on the production of leptoquarks as discussed in the following sections. 
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events 


selection steps (cumulative cuts) 
(NQl) I (NQ2) I (NQ3) I (NQ4) 


signal efficiencies: 


scalar LQ, tulq = 106.0 GeV 

56.1% 

49.7% 

41.8% 

(35.5 ±4.3)% 

scalar LQ, tulq = 190.5 GeV 

66 .6% 

61.7% 

56.6% 

(56.5 ±6.3)% 

vector LQ, tulq = 106.0 GeV 

55.1% 

48.8% 

42.5% 

(32.1 ±3.9)% 

vector LQ, tulq = 190.5 GeV 

59.0% 

55.9% 

51.3% 

(51.3 ±5.8)% 

expectee 

background events: 

e’''e“ —>■ 2 fermions 

17936 

6887 

19.1 

0.2 ±0.1 

e’''e“ —>■ 4 fermions 

5721 

3886 

226 

20.5 ±6.6 

two-photon 

14517 

209 

5.8 

< 0.4 

charged-current DIS 

16.9 

12.1 

8.8 

6.0 ±1.9 

total expected background 

38192 

10994 

259 

26.7 ±8.6 

o 

^served events: 

data 

43094 

11303 

292 

25 


Table 3: Selection efficiencies in the z/q channel. The listed efficiencies are for 
^/s^ = 200 GeV. Also shown are the remaining number of data events and 
the expected number of background events after each selection cut. All event 
numbers are quoted after a loose event preselection in which more than 2 tracks, 
more than 0.1 GeV total energy in the hadron calorimeter, | cos 0(^1 < 0.99, 
and 0.48 < < 0.52 is required. The errors are the total systematic 

uncertainties including the Monte Garlo statistical error. The component which is 
not due to Monte Garlo statistics is correlated between the separate contributions 
to the background. The determination of the systematic errors is discussed in 
Section |^. 
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5 Systematic Uncertainties 


Systematic uncertainties arise from (1) the luminosity measurement, (2) the lim¬ 
ited size of the simulated event samples, (3) the modelling of the signal processes, 
and (4) the modelling of the Standard Model backgrounds. 

The integrated luminosity has been varied by its uncertainty of 0.2%. The 
resulting error is negligible. The statistical error of around 1% on the signal 
efficiencies is taken into account. 


In the simulation of signal events described in Section |^, the leptoquark is 
assumed to decay before hadronisation. However, if the leptoquark mass and 
couplings to fermions are small, the hadronisation process could begin before the 
leptoquark decay. This leads to differences in the event properties, notably in 
the number of charged particle tracks per event. Similar to [Q, the dependence 
of the efficiency on the fragmentation modelling has been studied with a special 
version of PYTHIA where the two cases have been implemented [^. The resulting 
uncertainty in the efficiency is 5% for eq states and 11% for z/q states. Also, the 


parameters of the cone jet hnder have been varied (the minimum jet energy by 
(15 ± 5) GeV and the cone half angle by (1.0 ± 0.35) in the signal Monte Carlo. 
This results in an uncertainty of 1% for eq hnal states and 1% - 4% for z/q states. 


To assess the uncertainties in the modelling of the background, the simulated 
distribution of each quantity which is used in the selection (except for the variable 
1 i^ 2/23 + 1-4 ■ hi 1/12 in the z/q channel) is reweighted in turn to match the data 
distribution after all preceding cuts in the selection. The selection (cuts and 
likelihood) is left unchanged. The difference between the efficiencies for weighted 
and unweighted events is taken to be the systematic error. In principle, the 
reweighting procedure could give a larger systematic error if a signal is present 
in the data. This has been checked by adding signal Monte Carlo, normalized to 
a cross-section corresponding to the limit of our sensitivity, to the background 
Monte Carlo and performing the reweighting procedure again. No change in the 
derived systematic error was observed. 


To check the modelling of the variable In 7 / 23 -1-1.4-In 1 / 12 , hadronic two-fermion 
events from the high-energy data have been used. These events have been divided 
along the plane perpendicular to the thrust axis, and each half reconstructed 
separately. The distributions for such half-events have been compared between 
data and simulation, resulting in a 6 % systematic error. 


The systematic errors for all selection variables are added in quadrature. The 
total systematic error on the background is 31% for eq hnal states and 8% for 
z/q states. The largest contributions are the modelling of the quantity —qe cos 9^ 
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(16%) in the eq case and yu and 1/23 ( 6 %) in the uq case. 


The systematic errors on the efficiency are taken into account in the limit using 
the procedure of Highland and Cousins [^]. The background level is systemati¬ 


cally shifted down by the total systematic error ensuring conservative limits. 


6 Results 


Since no signihcant excess is observed, limits on the couplings are calculated for 
both scalar and vector leptoquarks and for those values of the branching ratio to 
electron-quark hnal states that are allowed in the case of chiral couplings. The 
leptoquark cross-section calculated with the program ERATO-LQ is taken as input 
to obtain limits on the couplings. The limit calculations are performed according 
to the procedure of which takes into account the expected background, the 
expected mass distribution, the signal efficiencies, and the observed candidates. 
The upper limits at the 95% CL of the coupling A as a function of the mass 
mLQ are given in Figures ^ and for the scalar and vector leptoquark states 
and for different j3e values. For /3e = 0 no production in eq collisions is possible. 
The limits for the states Si/ 2 (q=—2/3) and Vi/ 2 (q=—1/3) for /3e —> 0 are valid 
if the assumption of chiral couplings is dropped. If a coupling A = is 

assumed, where the electromagnetic coupling constant Oem is taken at the mass 
of the leptoquark, aem(^LQ) ~ 1/128, the mass limits range from 183 GeV to 
202 GeV depending on the leptoquark state. 

In this paper, only the case of chiral couplings is considered. If this constraint 
is dropped, then the constraints on the branching ratio shown in Table |l| no 
longer apply. For a given state, the resulting mass limits do not depend strongly 
on /3e. 


7 Conclusions 


We have searched for singly-produced leptoquarks in electron-photon interactions 
at e+e“ centre-of-mass energies between 189 and 209 GeV using data collected 
with the OPAL detector at LEP. No evidence is found for the production of 
these particles. Therefore, limits are set on the coupling A for scalar and vector 
leptoquarks as a function of the mass for different branching fractions /9e into 
eq hnal states. The A limits can be directly interpreted also as limits on A' for 
squarks in R-parity violating SUSY models with the direct decay of the squark 
into Standard Model particles. 
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Figure 1: Schematic diagram of the s-channel production of a leptoquark in 
electron-photon scattering. The photon is radiated by one of the LEP beams, 
fluctuates into a hadronic object, and one of the quarks interacts with an electron 
from the other beam. 
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likelihood output m^Q (GeV) 

Figure 2: Electron-quark decay channel; (a) distribution of the ratio E^j^fs^e 
after cut (EQ2); (b)-(d) distributions of the likelihood input variables —geCos6*e, 
cos6*g, and In 1/23 after cut (EQ4); (e) distribution of the likelihood output £ after 
cut (EQ4); (f) mass distribution of the selected leptoquark candidates after the 
full selection. The points with error bars are the data and the full line represents 
the total Standard Model background normalised to the data luminosity. The 
dashed and dotted histograms show the distribution for the scalar state S 1/2 at 
= 200 GeV with scaled masses of xlq = 0.53 and xlq = 0.95, respectively. 
The normalisation of the leptoquark signals is arbitrary. In plots (a) and (e), the 
arrow points into the accepted region. 
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Figure 3: Neutrino-quark decay channel: |cos 6 '(^)| distribution after cut 

(NQl); (b) distribution of the variable after cut (NQ2); (c) distribution 

of the variable In 1/23 + 1-4 In 1/12 after cut (NQ3); (d) transverse mass distribution 
of the selected leptoquark candidates after the full selection. The points with error 
bars are the data and the full line represents the total Standard Model background 
normalised to the data luminosity. The dashed and dotted histograms show the 
distribution for the scalar state S 1/2 at = 200 GeV with scaled masses of 
a^LQ = 0.53 and xlq = 0.95, respectively. The normalisation of the leptoquark 
signals is arbitrary. In plots (a), (b), and (c), the arrow points into the accepted 
region. 
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Figure 4: Upper limits at 95% confidence level on the coupling constant A for 
single production of the scalar leptoquark states for different branching fractions 

Pe- 
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Figure 5: Upper limits at 95% confidence level on the coupling constant A for 
single production of the vector leptoquark states for different branching fractions 

Pe- 
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